Superfluid pairing in a mixture of a spin-polarized Fermi gas and a dipolar condensate 
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ticated) laboratory conditions (see |25l - l27l | and references 
therein). As is well known, due to the density fluctua- 
tion of Bose condensates, mixing bosons can induce an 
effective attraction between two fermions [28.] . This was 
the mechanism in Efremov and Viverit's proposal |2^ 
to achieve p-wave Cooper pairings in a single-component 
Fermi-Bose mixture. A generalization of the scheme to 
a two-dimensional (2D) mixture involving dipolar bosons 
was recently carried out by Dutta and Lewenstein [30| . 
with the goal of realizing a superfluid with px + ipy sym- 
metry whose excitations are non-Abelian anyons that are 
the building blocks for topological quantum computation 
and Nishida and Tan [33] sought to 
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Abstract 

We consider a mixture of a spin-polarized Fermi gas and a dipolar Bose-Einstein condensate in which s-wave scattering between 
fermions and the quasiparticles of the dipolar condensate can result in an effective attractive Fermi- Fermi interaction anisotropic 
in nature and tunable by the dipolar interaction. We show that such an interaction can significantly increase the prospect of 
realizing a superfiuid with a gap parameter characterized with a coherent superposition of all odd partial waves. We formulate, 
in the spirit of the Hartree-Fock-Bogoliubov mean-field approach, a theory which allows us to estimate the critical temperature 
when the anisotropic Fock potential is taken into consideration and to determine the system parameters that optimize the 
critical temperature at which such a superfiuid emerges before the system begins to phase separate. 



I. INTRODUCTION 

Superfluid pairing of fermions in the Zth partial wave 
depends on the two-body scattering amplitude fi. Ac- 
cording to Wigner's threshold laws /; scales as (fcp.) 
for typical ground state atoms, where hk^p is the Fermi 
momentum. This reduces the critical temperature for 
the p-wave {I — 1) Bardee-Cooper-Schrieffer (BCS) su- 
perfluid to such a low level that it is virtually inacces- 
sible to current technology, except in situations where 
Wigner's threshold laws are not respected. A case in 
point is tuning atoms near the p-wave Feshbach reso- 
nance 043 1 where p-wave scattering can be dramati- 
cally enhanced. This excited the hope of realizing p- 
wave superfluids via p-wave Feshbach resonances in the 
ultracold atomic physics community. But, due to an in- 
creased number of inelastic collisions near the resonance, 
"Feshbach molecules" are short lived [B-Qj ^"^^ such a 
prospect appears difficult to attain. Another example, 
first pointed out by You and Marinescu [13], is subject- 
ing a dipolar Fermi gas to a sufficiently strong DC elec- 
tric field, where p-wave scattering can also be enhanced. 
A subsequent detailed investigation by Baranov et al. 
p~H indicates that the order parameter in such a spin- 
polarized dipolar Fermi gas is unusual in the sense that 
it is the superposition of all odd partial waves. Recently 
there has been an upsurge of effort aimed at achieving un- 
usual Fermi pairings in both two-component and single 
component dipolar Fermi gases [12 - lf8l ] along the direc- 
tion of You and Marinescu [13], motivated largely by 
recent rapid experimental advancement in achieving ul- 
tracold dense dipolar gases in ^'^K -^''Rb polar molecules 
[iiSl and in Cr [12, spin-1 Rb [2|, and Dy [H atoms. 

In the present work, we focus on an alternative route 
towards the same goal, seeking to engineer unusual super- 
fluid pairings in mixtures of different degenerate quantum 
gases. Cold-atom systems, besides being capable of un- 
precedented controllability (interaction strength, dimen- 
sionality, etc.), enjoy the distinct advantage over tradi- 
tional solid-state systems of being easily mixed to form 
new states of quantum gases and liquids. Bose-Bose, 
Bose-Fermi, and Fermi-Fermi mixtures of different iso- 
topes or species are now readily accessible under (sophis- 



achieve the same goal by inducing p-wave and higher 
partial-wave resonances from s-wave Feshbach resonance 
between atoms from different Fermi gases in different di- 
mensions. In a recent article 3J|, we explored singlet 
and triplet superfluid competition in a mixture of two- 
component Fermi and one-component dipolar Bose gases. 

In this paper, we consider a three-dimensional (3D) 
mixture between (nondipolar) fermions of mass mp in 
ground state \a) and (dipolar) bosons of mass mg in 
ground state \b), with their dipoles oriented along the 
z direction established by an external (either electric or 
magnetic) field. In contrast to the phonon spectrum of 
a nondiploar condensate, which is isotropic, the phonon 
spectrum of a dipolar condensate is anisotropic [ssl . IsgI ] . 
Such anisotropic phonons will thus induce an anisotropic 
Fermi-Fermi interaction, thereby providing an alterna- 
tive mechanism for creating a superfluid (in the mixture 
between a spin polarized Fermi gas and a dipolar conden- 
sate) with an order parameter which is the superposition 
of all odd partial waves. 

The outline of our paper is as follows. In Sec. |lTl 
we develop, within the framework of the Hartree-Fock- 
Bogoliubov mean-field approach, a theoretical foundation 
for self-consistently addressing challenges we will face in 
later sections. These challenges include the anisotropy of 
the induced interaction and the limitations due to phase 
separation. In Sec. IIIIl we calculate the critical tem- 
perature taking into consideration the anisotropic Fock 
potential with the help of the variational method, and 
we show that properly tuning the dipolar interaction can 
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dramatically increase the critical temperature in the long 
wavelength limit. In Sec. llVi we review the subject of 
phase separation and develop a procedure that allows 
us to determine the system parameters for achieving the 
optimal critical temperature before the system begins to 
phase separate. We summarize and conclude in Sec. |Vl 

II. THEORETICAL FORMULATION 
A. Bare Model 

In our model the low temperature physics arises from 
the interplay between short- and long-range two-body in- 
teractions. The short-range interaction is dominated 
by s-wave collisions which, given that the Pauli exclu- 
sion principle forbids such interactions taking place be- 
tween two identical fermions, are parameterized with 
two scattering lengths, ubb and ubf, characterizing s- 
wave scattering between two bosons and between a bo- 
son and a fermion, respectively. The long-range inter- 
action is the dipole-dipole interaction (restricted only to 
bosons in our model), which varies with r, the separa- 
tion between two particles in position space, according 
to C/oo (r) = (l — 3z^/r^) /r^, where (f represents 
the dipolar interaction strength. Assuming that the 
mixture is sufficiently large so that it can be treated as a 
uniform system with an effective volume V, we model 
such a mixture with a (grand canonical) Hamiltonian 
H = Hg + Hp + Hbf in momentum space according 
to 

Hb=J2 - ub) bibk + {2V)-^ X 

k 

J2 [UBB + UoDmbl^J^l^J^k'k, (1) 
k,k',q 

Hp =^{£k,F -up)aldw, (2) 

k,(T 

Hbf = UbfV^^ ^ a[ak+q^L^k'-q, (3) 

k,k',q 

where 6k (ok) is the field operator for annihilating a bo- 
son (a fermion) with kinetic energy ek,B = h^k^ /2mB 
(ek,F = fi.^fc^/2mi?) and chemical potential ub (mf); 
Ubb = Airh^aBB /rriB and Ubf — 4,TTh'^aBF /mBF with 
niBF = 2mBmp /{niB + mp) measure the related s- 
wave scattering strengths; Udd (k) = %-k(PP2 (cos6'k) /3 
is the dipole-dipole interaction in momentum space, with 
P2 {x) = {ix^ — 1) /2 the second-order Legendre polyno- 
mial and 0k (0k) the polar (azimuthal) angle of vector 
k. 



B. Effective (or Dressed) Model 

We consider the low temperature regime where vir- 
tually all bosons are condensed to the zero-momentum 



mode so that we can apply the usual Bogoliubov ansatz, 
in which 6k=o is treated as a c- number 6k=o- This trans- 
forms Eq. ^ into a diagonalized form with the famil- 
iar interpretation that the bosonic system described by 
Eq. Il]) is equivalent to a homogeneous dipolar con- 
densate with density ub = |6k=o| plus a collection of 
phonon modes that obey the Bogoliubov dispersion rela- 
tion 

Ek = VB hk ^l + (esfc)V2e<jd-P2(cos0k), (4) 

where Edd — 4:TT(P /{3Ubb) [131 measures the dipo- 
lar interaction relative to the s-wave collision, vb — 
yJuBUBB/'mB is the phonon speed, and = 
h/ y/^msnBUBB is the healing length. As can be seen, 
when Sdd > 1, phonons with fc — >■ acquire imaginary 
frequencies. In our study, then, we will limit Sdd to 
be less than one so that the dipolar condensate is stable 
against collapse. Integrating out the phonon degrees of 
freedom as in Ref. [28] leads to an effective Bose-Fermi 
model described by the Hamiltonian 

H^ = H"s+Hl (5) 

where 

H''b = V Qf/ssnl - UBUB^ (6) 
accounts for the condensate energy, and 

Hp = ^ (ek,_F - Up + UBpnB) a]^a\^+ 

k 

^ t^(q)«k+q4'-q«k'ak, (7) 

k,k',q 

is the effective Hamiltonian for the dressed fermions. In 
Eq. 0, 

f/(k)^ ^ (8) 

l + (^i3fc)' + 2eddP2(cos0k) 

is the phonon-induced Fermi-Fermi interaction in the 
static limit [2^ . It is well known that the dipolar in- 
teraction in position space is zero when averaged over all 
directions or, equivalently, the dipolar interaction in mo- 
mentum space, Udd (k), satisfies t/uu (k = 0) = and 
thus Udd (k) is non analytic when k approaches zero 
[ssf . For the same reason, the induced interaction U (k) 
is also non analytic when k approaches zero, but instead 
of being zero, J7 (k = 0) approaches the finite value 

U{Q) = -UIp/Ubb. (9) 

which is certainly different from [/ (k — >■ 0) according to 
Eq. due to its dependence on the polar angle in the 
denominator. As can be seen, the induced interaction 
depends on the dipole orientation differently than the 
direct dipole-dipole interaction and therefore provides an 
alternative mechanism for the exploration of the physics 
behind the anisotropy-related unusual Fermi pairings. 
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C. Hartree-Fock-Bogoliubov Mean-Field Theo- 
retical Treatment 



To study the superfluidity of a fermionic system de- 
scribed by the effective Hamiltonian ([7]), we use the 
self-consistent Hartree-Fock mean-field approach [39] in 
which we pair fermionic field operators in the two-body 
interaction part of the Hamiltonian in Eq. ([7]) according 
to Wick's theorem, associating direct pairing with the 
Hartree potential 

S]' = f/(0)nF, (10) 
exchange pairing with the Fock potential 

^(l^) = -^E^(k-k')(4'ak'), (11) 

and BCS pairing with the gap parameter between 
fermionic pairs of opposite momenta 



A(k) = i^[/^(k,k') 



ak'. 



(12) 



where np is the Fermi number density. In arriving at 
Eq. (I12p . we used a symmetrized induced interaction 



Ua (k, k') = i [C/ (k - k') - (7 (k + k')] , (13) 

in order to explicitly restrict the gap parameter to the 
sector in which A (k) = —A (— k), a constraint imposed 
by Fermi statistics (or equivalcntly the Pauli exclusion 
principle). This means that A (k) will be a superposition 
involving all odd partial waves when the Fermi-Fermi in- 
teraction is anisotropic fs^] , as is the case with our model 
here . 

The effective Hamiltonian in Eq. ^ now becomes the 
mean-field Hamiltonian 



H^' = Hi 



H 



F 1 



(14) 



where has already been defined in Eq. and the 
other two are defined, in terms of the self-consistent mean 
fields introduced in Eqs. ([T0| - (IT^ . as 



k 



[6 



S (k)] aj^Ok 



2^ 

k 



[A(k) aial, 



h.c 



(15) 



where is an effective chemical potential for fermions 
defined by 



fj,F — Up — [U (0) np + UBpnB] 



(16) 



and 



1 

k 



S](k)(4ak)-i5:A(k) (4at 



(17) 



or equivalcntly 



H''p = -V-U{Q)nl 



i^S(k)C/-i(k-k')S(k') 



k,k' 



l^A(k)C/^i(k,k')A* (k'). 



(18) 



k,k' 



where the inverse matrices are defined as 
F-^X^k" ^(k-k")?/"Mk" -k') = 4 k' and 
Ek" Ua (k,k") C/^i (k",k') = <5k,k'. 

Note that in problems where Bose and Fermi number 
densities are fixed, the second term in Eq. (jl6p (every- 
thing inside the square brackets) represents a constant 
shift and may thus be absorbed into the chemical po- 
tential Up, making it unnecessary to introduce ^p as we 
have done in Eq. p^ . This, however, no longer holds to 
be true in the case of phase separation. Phase separation 
redistributes the particle number densities among coex- 
isting phases so that the second term in Eq. and 
consequently ^p, will be different for different phases, 
even though the coexisting phases all share the same 
chemical potential up. Thus, in order for our theory 
to handle the subject of phase separation in a unified 
manner, we keep the second term explicitly. 

To obtain the partition function Z and the thermo- 
dynaniical potential $7, we make use of the Bogoliubov 
transformation between Ok and the fermionic quansi- 
particle operator bi^: ak = Ukfok + ^k&Lk and a^j^ = 
mJ^S^^ — Uk^k, where |uk|^ + l^kP = 1, and change the 
Hamiltonian Hp^ into the diagonal form 



2 ^ 

k 



[(a 



Sk + &L 



(19) 



where i?k — ■\/^J + A^(k) is the quansiparticle spec- 
trum of the BCS superfiuid. The corresponding grand- 
cannonical thermodynamic potential (density) [57 — 
- \nZ/{(3V) where 13 = 1/ {ksT) and ks the Bohzmann 
constant] is then given by 
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2V 



Ek) 



±J2^nf{-E^). (20) 



The phase transition from a normal gas to the BCS 
superfluid is known to be second order in nature. In 
the parameter regime near the second order phase tran- 
sition, the gap parameter is small. As a result, we may 
determine the properties of the Fermi gas close to the 
second-order phase transition by expanding the thermo- 
dynamical potential up to the fourth order in the gap 
parameter. We divide this expansion into two parts 



= f7o + ill. 



(21) 
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where 

^^0 = ^UBBn% - ubUb - (0) nj. 



i^ao(k) + i^^S](k){/-i(k-k')E(k'), 

(22) 



V 



k,k' 



represents the energy if the fcrmions were ah in the nor- 
mal state, and 

^1 = E ^ (1^) (k) SuM' + U^' (k, k')] A* (k') 



k,k' 



+ ^E"4(k)|A(k)|\ 



(23) 



represents the energy contributed by Cooper pairs, where 



ao (k) = ihi[/(-ek)] 



a2 (k) = 



1-2/ (a) _ tanhif^ 



26 



26 



"4 (k) = ^ 



1 - 2/ (6) 

2a 



/3/ (6) (1-/(6)) 



(24) 
(25) 

(26) 



are the relevant expansion coefHcients expressed in terms 
of the Fermi distribution function / (a;) = (e^^ + l) 
Equation (PT|) . (1^^ and (|23p serve as the foundation for 
our studies below. 



III. CRITICAL TEMPERATURE WHEN FOCK 
POTENTIAL IS INCLUDED 

In this section, we discuss how to compute the critical 
temperature when the Fock potential is included. The 
process begins with the equations for the Fermi density 



l^ifl-tanh 



V^2 



/?6 



and Fock potential 

^(k) = -^Et^(k-k')J 



1 — tanh 



/36 



(27) 



(28) 



which are obtained from np ~ —dfl/dfip and the saddle 
point condition for the Fock potential, 9r2/9E(k) = 0. 
The two equations are then coupled to the gap equation 



k' 

l^L/^(k,k')a4(k')|A(k')|'A(k'), (29) 



V 



due to the saddle point condition for the gap parameter, 
dfl/dA (k) = 0. In deriving Eq. (f29| . we renormalized 
the dipolar interaction by applying a procedure similar 
to the standard method for renormalizing the contact in- 
teraction ^] , but expressed in terms of vertex functions 
[Til, so that 



a2 (k) -^K{k) = 



tanh (/36/2) 
26 



Near the critical temperature, the gap is just about to 
vanish, and we can ignore the nonlinear term. This 
leads to the linearly coupled homogeneous equation 



A(k) = i^i?(k,k')A(k') 



(31) 



where i?(k,k') = ~Ua ik,k') K (k') is defined as the 
matrix element of R in momentum space. The critical 
temperature then corresponds to the largest temperature 
root of the characteristic (or secular) equation of the ma- 
trix R: 



Det (i? - /) = 0, 



(32) 



where / is the unit matrix. 

In general, one needs to include not only states on the 
Fermi surface but also those off the Fermi surface when 
solving simultaneously Eqs. ((27|) . p8| and (|32]) . which 
can make determining the critical temperature a compu- 
tationally expensive task. In order to gain quick insight 
into the critical temperature, we focus on the low tem- 
perature limit where we can make two significant sim- 
plifications. First, Eqs. (|27|) and (f28|) decouple from 
the gap equation, and second, mainly the states near the 
Fermi surface contribute to pairings. As such, the theory 
we formulate shall remain quantitatively accurate in the 
weak coupling regime and is expected to capture qual- 
itative features of the critical temperature under more 
general circumstances. 

To proceed, we first set the temperature in Eqs. (|27|) 
and (|28|) to zero. In this way, both the chemical and Fock 
potentials for a given n p can be obtained self-consistently 
from Eqs. p7| and (|28)) prior to being fed into the gap 
equation for estimating the critical temperature. Figure 
[1] (a) displays the chemical potential /i as a function of 
S (dots) when Sdd = 0.8 and A = 0.87 and Fig. [U (b) 
illustrates how the corresponding Fock potential E (k) 
changes with k and 0k when S = 0.8, where 



S = iBk% 



hj 7-1 



A^TTUBaBB 



and 



with 



X^N{e^^)Ulp/UBB = -^^^k% 

TT m^gp aBB 

s _ mpk'^ 



27r2ft2 



(33) 



(34) 



(35) 
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the density of states at the Fermi surface. The use of the 
superscript here is to stress that 

4 = (6^2^j.)'/^ , e% = h^ (4)' /2mF, (36) 

are, respectively, the Fermi energy and wavenumber of 
a spin-polarized non-interacting Fermi gas with number 
density np- Note that the rotational symmetry of the 
induced interaction with respect to the z axis implies 
that t/ (k — k') is a function of the relative azimuthal 
angle 0k — 0k' , and physical observables such as E (k) , 
being only dependent on the induced interaction after 
the azimuthal degree of freedom is integrated out in the 
manner 

/ C/(k-k')d(0k-0k'), (37) 
Jo 

are thus functions of (fc, 9^) only. As a result, 0k in S (k) 
[of Fig. lUb)] has been suppressed; this convention is to 
be extended to all the other observables. Additionally, 
for notational simplicity, C/ (k — k') and Ua (k, k') are to 
be understood, from now on, to represent the ones where 
azimuthal coordinates have been integrated out as in Eq. 
([57]). Figure 1(b) indicates that the Fock potential is 
lower in the external field direction (0k = and tt or 
the z-axis) than in the orthogonal direction ((?k = 7''/2) 
leading to a Fermi surface slightly elongated along the 
dipole orientation as shown in Fig. 1 (c) (dotted curve) , a 
feature apparently stemming from the anisotropic nature 
of the induced interaction in Eq. ([5]). 

To capture the anisotropic nature of the Fock poten- 
tial, we organize the single-particle dispersion 

in terms of two parameters, kp ((?k) and 7 (fc, 0^), which, 
in principle, can always be extracted from the numer- 
ically computed Fock potential. Physically, the for- 
mer defines the Fermi surface k = kp (dit), while the 
latter allows us to define mp^p (^k) as the renormalized 
mass along the 9k direction at the Fermi surface, where 
7f (^'k) = lir^k^kp (Sk) 7 (k, 6'k). These same physical in- 
terpretations also follow from the roles these quantities 
play in modifying the density of states. By definition, the 
number of states per unit volume, dN, along 6^ within 
the solid difi^erential angle dQ]^ — sin ^kC^^kC^^k, is given 
by 

- 7^ kT^^"'^ " ^ ^^^^ T^' (39) 
[2TTy |Vk4k| 47r 

where iV (ilk) is defined as the density of states. 
This means that in the limit where the Fermi 
surface varies slowly with 9^^ in the sense that 
\kp^ (6'k) dkp (6'k) /rf6'k| ^ 1, the density of states at the 
Fermi surface, 

TV (r!k) = N (eO ) [gp (9^) = ^p {9^) kp (9^) /4] , (40) 




FIG. 1: (a) The chemical potential, hf, is plotted as a func- 
tion of the unitless variable 5 = ^sfc^, where is the healing 
length [see Eq. H33p ]. The dotted curve is computed self- 
consistently from Eqs. (|27p and (|28p and may be compared 
to the solid one which is obtained variationally from Eqs. (|46p 
and (|54p . (b) The self-consistently computed Fock potential, 
S (fc,Sk), for 5 = 0.8. (c) A polar plot of the Fermi surface, 
kp (dk)- The dotted curve is extracted from the Fock poten- 
tial plotted in (b) and may be compared to the solid curve 
which is determined variationally from Eqs. (|46|) and (|51[) . 
(d) A polar plot of gp (^k), which is defined in Eq. (flO)) and 
is proportional to the density of states at the Fermi surface. 
The plot in (d) is determined variationally from Eqs. (|46p 
and (|53p . In (c) and (d), the polar angle is measured with 
respect to the positive fc^-axis. In each plot, we used A — 0.87 
and Edd = 0.8, where X — N (e^') Ugp/UsB is defined in Eq. 
(|34p . and Sdd = 'iTvd^ /{SUbb) measures the dipolar interac- 
tion relative to the s-wave interaction. 



simply corresponds to Eq. ([55)) for the isotropic case 
in which kp and mp are replaced, respectively, with 
kp (^k) and mpjp (9^). This result is in complete 
agreement with our expectation in light of the discussion 
above. Thus, we see that the Fock potential introduces 
anisotropy in the density of states so that the density of 
states at the Fermi surface is modulated along the angu- 
lar dimension by a factor of gp (^k)- 

In the low temperature limit ksT/e^p ^ 1, k'^K (k) is 
negligible except around the Fermi surface k = kp (^k), 
where it is sharply peaked compared to other momen- 
tum distributions inside the integration kernel in Eq. 
pip . This allows us to apply an analogous procedure, 
well-known in the study of isotropic BCS pairing [ioj . 
to single out the states near the Fermi surface as the 
dominant contribution to the gap equation. Indeed, 
we find, to a good approximation (correct up to a pre- 
exponential factor), that the gap at the Fermi surface 
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A {6k.) = A {k, 9k.) |fe=fcp(ek) obeys the equation 



A 61k In— 0— -^x 

47r Se^eT^ ^ 



where 
Ua{ 



9F {0k') Ua {0k, 0k') A (0k') sin0k'd0k', (41) 
Ua {k, 0k, k', Ok') \k=kF{0i^),k'=kF{ei^,) 



') 



U{0) 



(42) 

is the matrix element of the scaled interaction between 
two momenta on the Fermi surface. The critical tem- 
perature then becomes 



T 



Trks 



exp 



1 



(43) 



where lo is the largest positive eigenvalue of the eigen- 
value equation 



47r j gp {0k') Ua {0k, 0k') (^k') sin^k'd^k' = u^^ {0k) ■ 

(44) 

As can be seen, when we set gp {9k') = 1, we recover 
from Eqs. (|4T|) and (l44l) the corresponding equations in 
the absence of Fermi surface deformation. 

Inspired by the work in [4l|, [Z^] and the fact that our 
numerically determined Fermi surface fits well with an 
ellipsoid, we assume that the Femi particle number = 

a\jik/ distributes according to the variational ansatz 



nk-e 



(45) 



where Q [x\ is the Heaviside step function, and a is the so- 
called deformation parameter; the Fermi surface is pro- 
late when a < 1 and oblate when a > 1. An explicit 
evaluation of both kinetic energy and Fock potential in- 
dicates that a may be determined variationally by mini- 
mizing the (kinetic -f- Fock exchange) energy density: 



(46) 



where 



Y (a) = —a - 
^ ' 3 

4^ 



1 

3^ 



/ {9k, a, Edd, 5) sin 9kd9k- 



(47) 



InEq. dUl), 

I {9k,a,£dd,^) 



4d 

1 c 
16dd 




(48) 



plays the same role as the deformation function in dipolar 
Fermi gases [4ll | , but with a more complicated form that 
depends not only on the deformation parameter a but 
also on Sdd and 5, where 



c = 1 - Edd + 3e 
,52 



cos f k 
siv? 6'k + cos^ 9^ 



d = — {a^ sin^ 9^ + cos^ 61k) 



(49) 
(50) 



For an ellipsoid whose Fermi surface is defined by Eq. 
kp {9k) = k% {a-^ sin^ 6'k + cos^ 61k) ^^'^ 



(51) 

7_F (^k) = (a~^ sin" 6'k -f a' cos^ 6^kj , (52) 
and consequently according to Eq. (|40|) 

(53) 



gF(6'k) = (a"^sin2 6'k + a^cos2 6lk) . 



The corresponding chemical potential, calculated using 
fj-F ^ d£/dnp, is found to be 

MF 6 , . 1 /2 1 
4 =5^^")-5U" + ^ 
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A / ii:(6ik,a,edd,(5)sin6'krf6ik, (54) 
Jo 



where 



K {9k,a,edd,S) 
c 1 



= 3 



6 



c /'3c 24,, 



^3/2 
1 + 



tan 



(55) 



Figure 1(a) shows that Eq. ([Ml) produces a chem- 
ical potential (solid curve) that is in good agreement 
with the one obtained self-consistently from Eqs. ([27]) 
and ([25|) (dots). Figure 1(c) shows that compared to 
the self-consistent method (dotted curve) , the variational 
method (black curve) underestimates the prolateness of 
the Fermi surface, but nonetheless captures the essential 
anisotropic feature of the Fermi surface and does not, in 
our judgement, prevent us from using it as a computa- 
tionally economic tool for estimating the critical temper- 
ature. 

Our algorithm for computing the critical temperature 
is then as follows: we first calculate a by minimizing 
Eq. (1461) with respect to a; next we use it to fix 5^? {9^) 
according to Eq. ([5^: we then substitute gp {0k) into 
Eq. ([m to determine the eigenvalue w; finally we use 
Eq. ([^ to estimate the critical temperature. Figure 
[2| (a) compares the eigenvalue in the absence of dipo- 
lar interaction (dashed curve) with those in the presence 
of dipolar interaction (solid curves), indicating that in 
contrast to the former, which asymptotes to zero, in the 
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FIG. 2: The superfluid critical temperature, T, is propor- 
tional to exp(— l/uA) [Eq. (03])], where uj is the largest pos- 
itive eigenvalue of Eq. and A = N [e%) U%p/Ubb is 
defined in Eq. (|34l) . (a) The eigenvalue, uj, is plotted as a 
function of the unitless variable 5 — $,Bk%, where is the 
healing length [see Eq. H33p ]. The dashed curve corresponds 
to the absence of the dipolar interaction (edd = 0). The solid 
curves, from bottom to top, correspond to Sdd = 0.1, 0.2, 0.5 
and 0.8, where Edd = 4:Tvd^ /{3Ubb) measures the dipolar in- 
teraction relative to the s-wave interaction, (b) The critical 
temperature, T, in the units of Tp = ep/kB is plotted as a 
function of S. The curves, from bottom to top, correspond to 
£dd=0.6, 0.7 and 0.8. All curves are produced with A — 0.87. 



a much simplified f2, which, when expressed in terms of 
A= {Qt:'^Y^^ /'^''TTLp, takes the form 



1 o 



UBpnEnp 



5/3 



upnp, 



(56) 



where we have used Eqs. dH]), ^ and and 

ignored ili, the energy contribution from Cooper pairs, 
as we have done in other sections. 

Note that Eq. ([55)) has the same mathematical form 
as the one for a homogeneous Fermi-Bose mixture model, 
which has been extensively analyzed in [l^li^. Here, 
we recapture the main physics by performing the same 
analysis but in the chemical potential space [44]. Con- 
sider first the mixed phase where both ub and np have 
finite values. The corresponding chemical potentials are 
given by 



Ub 



Up 



duB 
dnp 



= Ubbtib + Ubftip, (57) 
= UBFnB + AnJ,^. (58) 



Eliminating the boson density, we arrive at a cubic equa- 



hniit of small (5, the latter approach finite values, which 
increase dramatically with the increase in Edd- This is 
explicitly shown in Fig. [2] (b), which is a plot of the crit- 
ical temperature for various Edd- Not shown in Fig. [Ifb) 
is the critical temperature in the absence of the dipolar 
interaction, which is orders of magnitude lower. 



IV. OPTIMAL CRITICAL TEMPERATURE 
PRIOR TO PHASE SEPARATION 

The issue of whether different phases can coexist in 
the same spatial volume (miscibility) or repel each other 
into separate spatial regions (phase separation) will in- 
evitably arise when one is considering systems made up 
of quantum gases with different species. The knowledge 
of phase separation is essential before one can apply the 
theory developed in the previous sections to estimate the 
critical temperature achievable before the system starts 
to phase separate. An understanding of the phase sepa- 
ration, even in view of the tremendous simplification due 
to the use of the variational method, would most likely 
remain a computationally expensive endeavor if one were 
to take into full consideration the Fock potential. Con- 
sequently, when studying phase separation we choose to 
ignore the momentum dependence of the induced inter- 
action by assuming U (k) — U{Q). The modification due 
to this assumption is not expected to undermine our goal 
(in this section) of gaining a qualitative understanding 
about the order of magnitude of the achievable critical 
temperature. This assumption allows the Fock poten- 
tial to cancel the Hartree potential identically, leading to 



tion for nj/^, 



Ubb 



+ Up 



Ubf 



UB = 0, (59) 



which allows us to determine Up given a set of up and 
Ub- Not all positive solutions to Eq. (15^1) are stable 

1 /3 

and therefore physically realizable; the values of rip are 
constrained by the stability condition which requires the 
following Hessian matrix to be positive semidefinite: 



\ drtpdns 

or equivalently 



dnsdrtF 
dn'i. 



1/3 



< A 



Ubb Ubf 
Ubf ^1^73 



2Ubb 

■i'-^BF 



(60) 



(61) 



There exists only one positive root to the cubic equation 
([59]) that satisfies the stability condition and this happens 
in the chemical potential space when 



Ubf 
Ubb 



Ub < Up < 



Ubf 
Ubb 



UB 



A'UU 



27 



^ BF 



(62) 



This means that there is a unique mixed state, which pre- 
cludes the possibility of phase separation involving more 
than one mixed phase. 

Next, in the region that supports a pure Bose phase, 
the effective chemical potential for fermions shall be less 
than zero, up — Ubp^ib < 0, which, when combined with 
Ub = Ubb^b for a pure Bose state, shows immediately 
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that the pure Bose phase exists in the chemical potential 
space with 



Up < Ubfub/Ubb- 



(63) 



Thus, we see that there is no overlap between the pure 
Bose and the mixed phase; the two phases are divided 
by the line up = UbfUb/Ubb, which is second order in 
nature. As a result, phase separation between the mixed 
state and a pure Bose state is impossible. 

The only remaining possibility is phase separation be- 
tween a mixed state with densities {npi, nBi) and a pure 
Fermi phase with densities {np2,nB2 = 0)- (A complete 
separation between fermions and bosons requires much 
higher densities than considered in the present work.) 
The mixed phase must share the same chemical and thcr- 
modynamical potentials with the pure phase, meaning 



UBpnBi + An 



2/3 _ ^„2/3 
2 . 5/3 _ 2 



Fl 



F2 ' 



(64) 



UBBriBi/'^ + UBpnBinFi + I^^fi' = ■^^'^F2^ (65) 



from which we find 



nBi 



U 



A 2/3/2 ^\ 



BF 



(66) 



1 /3 

where y = {nF2/'^Fi) is the sohition to the equation 



8y^ + 16y2 + 24y + 12 0, (67) 



where A is given by Eq. p4p except that e^;, is replaced 
with the Fermi energy e'pi of the mixed phase. The 
stability condition in Eq. (|6T|) then means that the stable 
mixture takes place when A < 1. 

Let us now turn to the question of how to identify the 
parameter space which optimizes the critical temperature 
before the system phase separates. A possible solution 
can be found in Ref. (29| . Here, we offer an answer from 
a different perspective. To be concrete, we consider a sys- 
tem in which tof = 6 amu, uib = 127 amu, ubb = 250ao 
(with amu being the atomic mass unit and oq the Bohr 
radius). Our task here is to find, for a given A and Edd, a 
set of {npUBTCLBF) which optimize the critical tempera- 
ture before the system phase separates. We begin with 
Eq. (j331), the equation for the critical temperature. We 
replace e^;. in Eq. (l43l) in favor of 5 defined in Eq. ([33)) 
and reorganize Eq. P5)) in terms of 5 as 



T 



-UBaBB 



kBiTip 



f {5) = 6^ exp 



1 



uj{S)X 



(68) 

where the argument in uj (6) is to stress that when A and 
Edd are fixed, the eigenvalue a; is a function of S only. 
On the one hand, in the limit of large S, w (6) is a small 
(positive) number [See Fig. [Sfa)], and thus T (6) will be 
low due to the exponential factor in Eq. (|68p . On the 
other hand, in the limit of small 6, T [5) is also low due 



to the factor 5"^ in Eq. ([68)1 . Thus, there will exist 
some 5 = Sp^ak at which T (6) reaches a maximum value, 
T (Speak)- Once S is fixed to Speak , the temperature T 
according to Eq. ([68]) increases with the boson density 
ub, and is highest when the largest ub is used. Pro- 
vided that Ub does not exceed values where three-body 
recombination begins to dominate the loss mechanism, 
the highest possible ub before phase separation equals 
UBi in Eq. ()66|) when npi = np. Finally, by combining 
Eq. dnni) with Eqs. dM]) and and solving them 

simultaneously, we find that the required {npUB^aBp) 
are given by 



np 



TLB 



\^ (mp/mBy 



^(^'^-%BSllakiy'-^f 

TT X^{mp/mB)'^ 

2(8^)4A(y^-i)' 



„ c2 ™-BF / 2 
ttBF = 2aBB0peak [V 

^ mp 



1 



(69) 

(70) 
(71) 



where y is the solution to Eq. (l67l) for the given A. The 
results of this optimization procedure are summarized in 
Fig. [3l A comparison between Figs, ^a.) and (b) indi- 
cates that mixing fermions with a dipolar condensate can 
result in a dramatic increase in the achievable tempera- 
ture T compared to mixing fermions with a nondipolar 
condensate; the order of magnitude in T can be 100 nK 
for the former and only a fraction of 1 nK for the lat- 
ter. Figures [3jc) and (d) display, respectively, the ubf 
and (nB.Tip) required to realize the temperature in Fig. 
[3l[b) when £dd = 0.8 (solid curve). For example, when 
A = 0.874, a system with obf = 275ao, ub = 2.76 x lO^'' 
cm"'^ and np = 6.88 x 10"'^^ cm~^ has a critical temper- 
ature of 80 nK. 



V. CONCLUSION 

The ability to easily mix cold atoms of different species 
opens up another exciting avenue for engineering quan- 
tum systems with novel properties. Mixing a dipolar 
condensate can induce, between two fermions in a Fermi 
quantum gas, an effective attractive interaction which, 
owning to its anisotropic nature, can dramatically en- 
hance the scattering of fermions in the odd-parity chan- 
nels. We have explored such an interaction for the pur- 
pose of achieving, in a mixture between a spin-polarized 
Fermi gas and a dipolar condensate, a superfiuid with a 
gap parameter characterized with a coherent superposi- 
tion of all odd partial waves. We have focused our effort 
on determining the critical temperature, a task far more 
challenging for our model than for pure fermionic systems 
since, first, adding dipolar bosons to fermions dramati- 
cally increases the system parameter space and, second, 
the system may phase separate. We have formulated, 
in the spirit of the Hartree-Fock-Bogoliubov mean-field 
approach, a theory which allows us to treat calculating 
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the critical temperature and phase separation in a unified 
manner. We have apphed this theory to estimate the 
critical temperature when the anisotropic Fock potential 
is taken into consideration and to identify the parameter 
space which optimizes the critical temperature before the 
system begins to phase separate. 
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FIG. 3: (a) In the absence of the dipolar interaction {edd = 0), 
the superfiuid critical temperature, T, is plotted as a function 
of the unitless variable A, where X = N (e'^) Ugp/UsB is de- 
fined in Eq. (|34|) . (b) For a nonzero dipolar interaction the 
critical temperature, T, is plotted for Sdd = 0.6 (long dash), 
0.7 (short dash), and 0.8 (solid), where Sdd = 47rd^/(3?7ss) 
measures the dipolar interaction relative to the s-wave inter- 
action, (c) aBF and (d) densities ub [black curve in (d)] and 
np [gray curve in (d)] that are required to produce the (op- 
timal) critical temperature for Sdd ~ 0.8, represented by the 
solid curve in (b). The black dots mark the points at which 
ub has reached 5 x 10^* cm~^. All curves are produced with 
rriF = 6 amu, niB = 127 amu, and ubb = 250ao. 
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